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ABSTRACT 


The  problem  of  radiation  from  sources  in  the  presence  of 


smooth,  convex 


new  techniques  based  on  the  spectral  domain  approach  and  an 


current,  also  are  discussed.  Some  numerical  results  derived  from 


the  spectral  domain  formulas  are  presented  and  a comparison  with 


available  theoretical  and  experimental  data  is  included 
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1.  Introduction 


The  problems  of  radiation  from  sources  in  the  presence  of  impenetrable 
smooth  convex  objects  and  the  diffraction  of  a plane  wave  by  such  objects  are 
of  great  practical  interest  in  the  design  of  antennas  on  structures,  e.g., 
conformal  arrays.  Unfortunately,  the  exact  analytical  solutions  to  these 
problems,  based  on  the  methods  of  "separation  of  variables"  or  "function- 
theoretic"  procedures  (Wiener-Hopf  technique,  residue  calculus,  etc.),  exist 
only  for  a very  limited  number  of  scattering  geometries.  Furthermore, 
the  exact  solutions  are  typically  highly  complex  in  nature;  hence, the  process 
of  extracting  numerical  results  from  them  can  be  very  time-consuming  and  is 
by  no  means  trivial.  This  situation  has  motivated  many  researchers  to  explore 
approaches  to  the  problems  of  radiation  and  scattering  from  smooth  convex 
structures . 

In  the  low  and  resonant  frequency  ranges,  several  reliable  numerical 
procedures,  e.g.,  the  moment  method,  are  available  for  solving  the  radiation 
and  scattering  problems.  However,  in  the  high  frequency  domain,  numerical 
techniques  based  on  matrix  methods  become  unwieldly  if  not  impractical,  prompting 
one  to  employ  asymptotic  techniques  suitable  for  large  k(»27r/X),  where  X is 
the  wavelength  of  the  illuminating  wave. 

In  this  work,  we  begin  by  presenting,  in  Sec.  2 , a survey  of  various 
high  frequency  asymptotic  techniques  for  the  problem  stated  above.  The  survey 
will  be  necessarily  brief,  and  will  cover  only  the  highlights  of  a number  of 
important  approaches  to  the  problem  at  hand,  viz.,  Fock's  theory,  the  geometrical 
theory  of  diffraction  (GTD) , and  the  direct  integral  equation  approach.  The 
reader  interested  in  further  details  may  choose  to  consult  the  works  of  Bowman, 
et  al.  [1],  Uslenghi  [2],  and  Kouyoumjian  [3]. 
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In  Sec.  3,  we  consider  the  generalization  of  GTD  and  present  some  new 


approaches  to  the  curved  surface  radiation  and  scattering  problems.  Some 
numerical  results  based  on  one  of  these  new  approaches  are  presented  in  Sec.  3 
and  a comparison  with  other  available  methods  are  included. 

2.  Survey  of  Available  High-Frequency  Asymptotic  Techniques 
2.1  Watson  Transformation 

One  of  the  first  successful  attempts  to  derive  an  asymptotic  expansion 
for  the  far-field  generated  by  a point  source  located  in  the  proximity  of  a 
conducting  surface  was  made  by  G.  N.  Watson  in  1918  [4].  His  method, 
essentially,  consisted  of  two  steps:  1)  transforming  the  original  infinite 
series  solution  into  a contour  integral  (by  Cauchy's  residue  theorem);  2) 
deforming  the  contour  of  integration  so  as  to  capture  a set  of  complex  poles 
of  the  integrand.  The  original  integral  is  then  expressed  in  terms  of  an 
infinite  series  which  converges  very  rapidly,  provided  the  observation  point 
is  in  the  shadow  region.  The  first  few  terms  of  this  series  were  later  inter- 
preted as  "creeping  waves."  The  method  was  first  applied  to  a sphere  and 
circular  cylinder,  and  later  to  some  other  geometries  as  well.  The  mathematical 
rigor  of  the  method  was  the  subject  of  further  investigations  by  other  researchers 
([5],  [6]  and  [7]).  Although  Watson  transformation  can  only  be  applied  to  a 
few  simple  goemetries,  e.g.,  the  sphere,  cylinder,  cone,  spheroid,  it  is 

still  regarded  as  one  of  the  cornerstones  of  the  more  general  high  frequency 
techniques  because  of  its  mathematical  rigor.  Watson  transformation  is 
especially  powerful  in  the  shadow  region  of  the  geometric  optics  field.  In 
the  lit  region,  the  above-mentioned  contour  integral  is  evaluated  using  the 
"stationary  phase"  method  and  yields  the  reflected  field  from  the  surface. 

In  this  region,  the  most  significant  contribution  to  the  total  scattered  field 
typically  comes  from  the  surface  current  induced  on  the  smooth  convex  part  of 
the  object;  the  so-called  "Physical  Optics"  approximation  can  be  applied 
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([8],  [9],  and  [10])  to  derive  the  reflected  field.  The  Physical  Optics 
method  is  based  upon  approximating  the  induced  surface  current  in  the  lit 
region  of  the  object  by  the  current  that  would  be  induced  on  the  local  tangent 
plane,  and  by  assuming  that  the  surface  current  is  zero  in  the  shadow  region. 

The  far  field  is  constructed  by  substituting  the  above  estimate  for  the 
induced  surface  current  in  the  integral  representation  of  the  scattered  field, 
and  evaluating  the  same  in  an  asymptotic  sense.  The  dominant  term  of  the 
asymptotic  expansion  of  this  integral  can  be  shown  to  be  identical  to  the  first  term 
of  the  Luneberg-Kline  expansion  of  the  geometrical  optics  far  field  ([11]  and 
[12]).  However,  the  higher-order  terms  derived  from  the  physical  optics 
approach  do  not  provide  us  with  the  correct  result  in  the  shadow  or  transition 
regions  where  the  diffracted  field  contributes  the  most. 

In  the  next  section,  we  discuss  Fock's  theory,  which  can  fill  the  gap 
between  the  Physical  Optics  in  the  lit  region  and  the  "creeping  wave" 
representation  in  the  shadow  region. 

2.2  Fock's  Theory 

The  region  between  the  lit  and  the  shadow  part  on  a surface  is  called 

"penumbra  region."  The  angular  width  of  this  region  is  approximately  given 
1/3 

by  (Xr^/ir)  where  X is  the  wavelength  of  the  illumination  and  rQ  is  the 
radius  of  curvature  of  the  surface  of  the  object  in  this  region  in  the  incident 
plane  (Fig.  1).  Fock's  theory  invokes  the  principle  of  loaal  character  of 
the  field  in  the  penumbra  region  [13]  and  is  based  on  the  conjecture  that  all 
bodies  with  a smoothly  varying  curvature  have  the  same  current  distribution 
in  the  penumbra  region,  provided  that  the  curvature  and  the  incident  wave  are 
the  same  near  the  point  under  consideration.  This  principle  allows  one  to 
locally  replace  the  surface  of  the  object  by  a portion  of  a paraboloid  of 
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Figure  1:  Section  of  the  body  in  the  plane  of  incidence. 


Shadow 

boundary 


Figure  2:  Plane  wave  incident  upon  a smooth  convex  body. 


revolution.  A unique  feature  of  the  expressions  for  Fock  currents  is  that 
they  provide  a convenient  transformation  of  the  geometric  optics  currents 
in  the  lit  region  into  the  creeping  wave  currents  in  the  shadow  region. 

Fock  himself  deduced  the  pertinent  formulas  for  the  surface  currents  by 
treating  a convex  body  problem  [14]  described  below. 

Consider  a convex  body  and  a plane  wave  incident  in  the  direction  of  the 
x-axis.  If  the  equation  of  the  surface  is 

f (x,y,z)  - 0 (1) 

then  the  curve  representing  the  boundary  of  geometrical  shadow  is  given  by 

f(r)  - 0 , ||  » 0 (2) 

Consider  a point  0 on  the  boundary  of  a shadow  region  where  we  set  up 
a rectangular  coordinate  system  as  shown  in  Fig.  2 (z:  normal  to  the  surface, 
x:  in  the  direction  of  propagation,  and  y is  the  tangent  to  the  boundary  of 
shadow).  In  the  vicinity  of  this  point,  the  surface  of  the  body  could  be 
locally  replaced  by  a paraboloid  of  revolution  which  is  expressed  by  the 
equation. 

z + 1/2  (ax2  + 2bxy  + cy2)  ■ 0 (3) 

Each  of  the  field  components  satisfies  the  Helmholtz  equation 
(V2  + k2) y =*  0 (4) 
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The  fact  that  the  incident  wave  travels  along  the  x-axis,  suggests  that  ¥ 


be  written  in  the  form 

S'  =■  fe"jkx  (5) 

where  an  exp (jut)  time  dependence  has  been  assumed.  Substituting  (5)  in  (4) 
gives 

V2Y  - 2jk  4~  V =■  0 (6) 

dX 


i 


At  this  point,  two  basic  assumptions  are  introduced  in  Fock's  theory, 

viz. 

i)  f's  are  relatively  slowly  varying  function  of  coordinates 

ii)  S'  varies  more  rapidly  in  the  z-direction  than  in  x and  y, 

i.e. , 


3^ 

3z 


o(£  5), 

m 


3Y 

3x 


o(—  ?)  il 

3y 


0 (-  S') 

m 


Based  upon  (7),  we  can  write  (6)  as 


0 


(7) 


(8) 


and  consequently  m'  * m (m  is  very  large),  where  the  terms  of  order  1/m 
have  been  omitted. 

Inserting  these  estimates  and  assumptions  into  the  Maxwell's  equation, 

we  can  find  some  simple  expressions  for  all  the  field  components  in  terms  of 

H and  H . If  we  write  H as 
y z y 
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H - H°  e"jkx  f 

y y 


where  is  the  magnitude  of  the  incident  wave  at  infinity,  then  V must 
satisfy 


CIO) 


with  boundary  condition 


£ U/  

— -jk  (ax  + by  ) ¥ 


1 

% 


(ID 


on  the  surface  of  the  body.  Eqn.  (11)  is  the  simplified  version  of  the 
Leontovich  boundary  condition  where 


n * e - j 


4ttct 

ek 


The  final  solution  for  on  the  surface  of  the  body  which  satisfies 
the  boundary  condition  and  the  condition  at  infinity,  may  be  written  in  the 
form 

= H®X  G(£,q)  (12) 

sx 

where  H » external  field 


G(£,q)  = e“(j/3)  S V^C.q) 


V^(C,q)  ■ Fock  function  defined  in  the  Appendix  A. 

5 - m(ax  + by)  ■ reduced  distance  from  the  shadow  boundary  = i/d. 

2 1/3 

d ■ the  width  of  penumbra  region  ■ (2r  / k) 

o 

l m distance  between  the  observation  point  and  the  shadow  boundary  along  the 
incident  ray  (Fig.  3). 
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Figure  4:  Coordinates  of  observation  ooint  in  tertis  of  l 


q - -jin/ /n*  -(j/n)  \h~  (■  0 for  conducting  body) 


The  other  tangential  component  of  the  magnetic  field  11^  on  the  surface 
of  the  body  can  be  obtained  in  a similar  manner 


H *■  HeX  [-  j-  e”(j/3)C  /3  f(C)] 
x z m 


(13) 


where  f(C)  is  another  Fock  function  defined  in  Appendix  A.  Fock's 
formulas  not  only  give  the  surface  value  of  the  field,  but  also  can  be 
utilized  to  find  the  field  in  the  proximity  of  the  object.  For  a plane  wave 
incidence,  the  first  order,  i.e.,  0(l/m)  terms  for  the  scattered  field  within 
a certain  layer  around  the  object,  can  be  written  as 


H - 0,  H = H°  e"jkX  * (C,C),  H = H°  e~jkx  $(£,n) 
x y y z z 


E -(j/m)H°  e~iVx  3Y/3C,  E„  * H_,  E_  * -H 
x y 


(14) 


y z z y 


where 


2,2  2 

5 - 2am  [z  +(l/2)(ax  + 2bxy  + cy  ) ] = reduced  height  from  the  surface 


of  the  body  (Fig-  4) 


¥ = -je 


(JCC-J/3  C3)  r -jCt  , wJ_(t)-qWl(t) 


/ twl(t-0  - ^ 


2rt)-qw2(ty  W2(t-^)]  dC 


(15) 


- <1/3)5*  ,ct  w (t) 

/ «j5 


The  path  of  integration  for  $ and  V is  shown  in  Fig.  5 . 


dt 
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Fock's  expressions  for  the  field  components  in  the  penumbra  region 
(£  - 0)  can  be  extended  to  the  shadow  region,  by  introducing  some  modifications 


in  the  definition  of  parameter  £. 


Goodrich  [15]  has  generalized  the 


argument  used  by  Fock  in  the  penumbra  region  to  anywhere  in  the  shadow 
region  by  introducing  a new  set  of  variables,?  and  £,for  the  incremental 
distances  along  the  path  leading  into  the  shadow  region.  In  this  generalization, 
the  parameter  £ as  defined  in  (12)  is  replaced  by 


f /kR(s) . 1/ 3 ds 

i K 2 > R(s) 


where  s is  the  arc  length  along  the  geodesics  which  originate  from  the  shadow 
boundary  and  go  into  the  shadow  region  along  the  surface,  and  R(s)  is  the 
radius  of  curvature  of  the  surface  along  the  geodesics.  For  the  case  of  a 
circular  cylinder  of  radius  a (Fig.  6),  the  expression  £ simplifies  to 


£ - (ka/2)1/3  9 - s/d 


Fock  also  treated  the  case  where  the  point  source  was  very  close  to 
the  surface  of  the  body.  He  analyzed  the  radiation  of  electric  dipoles  near  a 
spherical  model  of  the  earth  [16]  and  derived  the  formulas  for  the  scattered 
fields  in  terms  of  functions  (attenuation  functions)  similar  to  4)  and 
V , which  are  valid  both  in  the  shadow  and  transition  regions  [17].  Fock's 
assumptions  were  later  proven  in  a more  systematic  and  mathematically  rigorous 
manner  by  Cullen  [18]  and  Hong  [19]  by  using  a direct  integral  equation 
approach.  This  method  is  described  in  the  next  section. 
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where  n(r)is  the  outward  unit  normal  to  the  surface  at  r,  H^"nc(r)  is 
the  incident  magnetic  field  on  the  surface  (S)  of  the  body, and  R ■ r - r 
(r ' is  a variable  point  on  the  surface) . 

Fock  used  this  integral  equation  to  deduce  the  important  principle 
of  local  character  of  the  field  in  the  penumbra  region.  Cullen  derived  a 
first-order  asymptotic  solution  to  (18)  which  agreed  with  Fock's  results 
given  in  (12)  and  (13).  Cullen's  method  consists  of  transferring  the  two- 
dimensional  integral  equation  (18) , in  the  penumbra  region,  to  a one-dimensional, 
Volterra-type  equation.  This  is  accomplished  by  applying  the  stationary 
phase  technique  to  the  original  integral  while  integrating  with  respect  to 
one  of  the  variables.  The  resulting  one-dimensional  Volterra  equation  is  then 
solved  in  Cullen's  method  by  the  Fourier  transform  technique.  A similar 
procedure  was  used  by  Hong  [19]  to  analyze,  asymptotically,  the  diffraction 
of  electromagnetic  and  acoustic  plane  waves  by  smooth  convex  bodies.  We  will 
now  proceed  to  explain  Hong's  method  in  a little  more  detail  by  referring 
back,  once  again,  to  the  integral  equation  (18).  The  surface  is  parametrized  by 
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the  geodesic  coordinate  system  (o,v)  such  that  the  shadow  boundary  for  the 
incident  plane  wave  traveling  along  the  tangent  a(0,v)  to  the  v = 0 curve  is  the 
a ■ 0 curve.  The  quantities  a(a,v),  b(a,v)  and  n(o,v)  form  a right-hand  local 
orthonormal  basis  (n  = a * b)  (Fig.  7). 

Since  the  incident  field  has  a phase  factor  e r(af0)^ 

we  write  the  surface  current  in  the  form 


J(r)  - [Ia(r)  a (r)  + Ife(r)  b(r)]  e 


-jka 


(19) 


where  a is  the  arc  length  along  the  geodesic.  Substituting  (19)  back  into 
(18)  and  restricting  the  resulting  equation  to  the  points  on  the  geodesic  v»0, 
we  obtain  two  coupled,  two-dimensional  integral  equations  for  l^(a,0)  and 
I (a,0).  It  can  be  shown  that  these  integrals  have  saddle  points  at  v*0  (for 

D 

the  v- integration) . Applying  the  "steepest  descent  path"  method  to  v-integration, 

2 1/3 

and  keeping  the  terms  up  to  the  order  l/M^,  where  = (kp^  (o,0))  , we 

obtain  the  following  decoupled  one-dimensional, Volterra-type  integral  equations 

for  I,  (5,0)  and  1(5,0) 
b a 


Ia(C,0)  - 2 ljnc(£,0)  - / dt  Ia(S,C)  Kct($-t)  +0(MqJ) 


“3. 


Ib(5,0)  - 2 ljnc  (5,0)  - / dx  Ib(C,0)  K^S-t)  + 0(M“J) 


,~3, 


(20) 


p (cr,v)  is  the  radius  of  curvature  of  the  surface  along  geodesics  (v  = constant 
o 

curves)  at  point  (a,v). 

Solving  (20)  by  Fourier  transforms,  we  obtain  the  expression  for  the 
induced  currents  in  the  penumbra  and  shadow  regions, and  the  first-order  solutions 
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are  found  to  be  the  same  as  those  of  Fock  and  GTD  [20].  One  of  the  important 
conclusions  drawn  from  Hong's  solution  is  that  the  leading  term  in  the 
asymptotic  expansion,  which  is  the  same  for  the  acoustic  and  electromagnetic 

I problems,  is  independent  of  curvature  in  the  direction  transverse  to  the 

geodesic,  provided  the  divergence  factor  is  suppressed.  However,  we  should 

bear  in  mind  that  Hong's  method  was  designed  for  the  case  of  axial  incidence 

on  symmetric  objects, and  in  this  case,  the  geodesics  are  torsionless.  The 

above  conclusion  does  not  seem  to  be  valid  in  the  cases  where  the  rays  have 

nonzero  torsion  ([21],  [22]).  In  Hong's  expressions  for  the  surface  current, 

the  transverse  curvature  has  only  a second-order  effect.  It  was  also  shown 

-2/3 

that  up  to  the  terms  of  order  (kp^)  in  the  asymptotic  expansion,  the 
tangential  and  binormal  components  of  the  creeping  waves  are  not  coupled. 

Both  Fockls  theory  and  the  "direct  integral  equation  approach"  give  the 
induced  surface  current,  or  the  scattered  field  in  the  neighborhood  of  the 
surface  of  the  scatterer,  due  to  an  incident  plane  wave.  These  expressions  can 
also  be  used  to  derive  the  radiated  field  via  the  use  of  the  reciprocity  theorem. 

The  methods  which  have  been  discussed  thus  far  are  mathematically 
rigorous.  However,  they  are  limited  in  the  scope  of  their  application  to 
geometries  satisfying  some  special  smoothness  and  symmetry  criteria.  "Geometrical 
theory  of  diffraction"  (GTD),  which  we  discuss  in  the  next  section,  has  a 
broader  scope,  although  it  does  lack  the  mathematical  rigor  of  approaches 
described  until  now. 

2.4  Geometrical  Theory  of  Diffraction  (GTD) 

Geometrical  theory  of  diffraction  (GTD),  developed  by  J.  B.  Keller 
([20],  [23],  [24],  [25],  and  [26]),  is  a generalization  of  geometrical  optics. 
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It  is  based  upon  the  assumption  that  fields  propagate  along  rays.  Keller's 
major  contribution  was  to  introduce  the  new  kinds  of  rays  called  the  "diffracted 
rays"  which  together  with  the  geometrical  optics  rays, constitute  the  total 
field.  In  our  problem,  viz.,  source  radiating  in  the  proximity  of  the  smooth 
object,  the  diffracted  rays  travel  along  the  curves  on  the  surface  of  the 
scatterer.  By  applying  Fermat's  principle  to  these  surface  rays,  we  conclude 
that  the  above-mentioned  curves  should  be  geodesics  on  the  surface  of  the  body. 

In  the  GTD  procedure,  one  assigns  a value  to  the  field  along  each  ray  of  these 
surfaces.  The  total  field  at  any  point  in  the  space  is  the  sum  of  the  fields 
due  to  various  rays  (incident,  reflected  and  diffracted)  passing  through  that 
point.  An  important  advantage  of  the  GTD  approach  is  that  it  can  be  applied 
to  both  scalar  (acoustic)  and  vector  (electromagnetic)  problems  and  to  smooth 
convex  objects  of  an  arbitrary  shape. 

Consider  the  problem  of  determining  the  radiated  field  of  a scalar  point 
source  located  on  the  surface  of  a smooth  convex  opaque  body.  If  the  observation 
point  is  in  the  shadow  region,  the  ray  paths  originating  at  Q and  reaching  P 
(observation  point)  are  comprised  of  two  sections.  One  of  these  sections  follows  the 
straight  line  path  P^P,  while  the  other  travels  along  a geodesic  on  the  surface 
(Fig.  8).  Let  us  consider  the  propagation  of  the  field  along  each  section  separately. 

a)  Rays  in  free  space:  Behavior  of  the  fields  along  these  rays  can  be 
determined  by  obtaining  a high-frequency  asymptotic  solution  to  Maxwell's 
equation  in  a source-free  homogeneous  isotropic  medium.  We  begin  with  the 
Luneberg-Kline  asymptotic  expansion  of  the  electric  field  ([11]  and  [12]): 


and  insert  it  into  the  Maxwell's  equations.  This  results  in  the  following 
equations  governing  the  propagation  of  electromagnetic  fields  along  the  rays. 


[VS(r)]  = 1 (Eikonal  equation) 


2(7S*7)e  + (V^S)e  = e , (Transport  equation)  (23) 

m m m-i 


VS'e  = -7-e  , (Gauss's  Law) 

m m— i 


e_^  =0,  m*0,  1,  2,  ... 


The  zeroth-order  solution  to  the  above  system  of  equations, which  turns  out 
to  be  in  agreement  with  what  one  would  obtain  by  geometric  optics,  may  be 
written  as 


E(o)  - EQ(0)e 


-jk0S(0) 


/ 

(Pj^+cO  (p2+a) 


where  a is  the  distance  traveled  along  the  ray  from  the  reference  point 
0(j=0)  on  the  ray  path,  and  p2  are  the  principal  radii  of  curvature  of 

the  wavefront  at  a=0.  It  is  apparent  that  the  expression  fails  when  a^-o^  or 
o = ~P2»  i.e.,  at  the  caustic  lines  (Fig.  9).  In  the  cases  where  it  is 
convenient  to  choose  the  point  of  diffraction  on  the  surface  of  a body  as 
the  reference  point  0,  the  formula  (2.5)  should  be  modified  as  follows 


o (p+o ) 
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In  these  cases,  the  point  of  diffraction  itself  is  a caustic,  and  p is  the 


distance  between  this  point  and  the  second  caustic. 

b)  Surface  Rays:  These  rays  follow  the  surface  S along  the  geodesics 
into  the  shadow  region,  and  shed  off  energy  tangentially  as  they  propagate. 

In  order  to  study  the  behavior  of  the  field  along  these  rays,  we  introduce  a 

* A A 

special  ray-fixed  coordinate  system, a, n,b. 

a:  Unit  vector  tangent  to  the  ray;  n:  outward  unit  normal  to  the 
surface;  and  b = t * n or  binormal  direction;  a vector  field  can  be  decomposed 
into  its  components  along  these  unit  vectors  as 

E - E^a  + E^n  + Ebb  (27) 

At  this  point,  several  important  assumptions  are  introduced  in  the  GTD 
approach  [20]: 

i)  E and  H are  orthogonal  to  each  other  and  to  the  ray. 

ii)  Variation  of  the  phase  of  the  field  along  the  ray  is 
the  same  for  both  fields. 

iii)  and  E^  propagate  independently,  and  E^  = 0. 

iv)  E^  satisfies  the  scalar  wave  equation  (V2+k2)  u=«0  with  the 

boundary  condition  u=*0  on  the  surface  S,  while  E^  satisfies 

the  same  equation  with  the  boundary  condition  ■ 0. 

on 

The  next  step  in  the  GTD  approach  is  to  conjecture, on  the  basis  of  the 
solution  to  some  canonical  problems, that  the  surface  field  propagating  along 
each  ray  is  comprised  of  an  infinite  set  of  "modes."  Along  a ray-fixed  path 
GTD  assigns  a complex  value  to  each  component  of  the  field  associated  with 
the  individual  modes.  The  propagation  of  these  modal  field  is  described  by 
the  equation 
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j (*0-ko) 


a(a)  ■ A (a)  e 


(28) 


when  a is  the  distance  between  an  arbitrary  point  along  the  ray  and  the  source 
Q and  is  the  phase  of  the  field  at  the  source  point.  Next,  invoking  the 
principle  of  conservation  of  energy  between  two  adjacent  rays,  and  using  the  fact 
that  the  surface  rays  shed  energy  off  tangentially,  we  can  arrive  at  the 
following  expression  for  a(a) 


a(a) 


exp  [ - j ka 


(o' )da ' ] 


(29) 


where  a (a)  is  the  "attenuation  constant,"  K is  proportional  to  the  strength  of 

1/2 

the  source,  and  , df^  and  p are  shown  in  Figure  10.  The  quantity  [dV^/  (pdV^)  ] 
indicates  the  spreading  of  the  surface  ray  tube"  as  it  travels  along  the  surface. 
Equations  (26)  and  (29)  describe  the  laws  of  propagation  for  the  rays  which 
originate  from  the  source  point  Q,  are  diffracted  at  P^,  and  reach  the 
observation  point  P.  To  complete  the  solution,  we  need  to  determine  the  actual 
values  of  the  fields  from  these  equations.  These  require  the  knowledge  of  <5^ 
and  K,  which,  in  turn,  are  related  to  the  initial  values  of  the  rays  QP^  and 
P^P  as  well  as  the  attenuation  constant  a (a).  The  initial  value  of  the  field 
at  0 is  related  to  the  strength  of  the  source  by  L(Q),  the  so-called  "launching 
coefficient,"  while  the  initial  value  of  the  field  at  P^  is  related  to  the 
actual  field  on  the  surface  at  P^  through  the  "diffraction  coefficient"  D(P^). 

If  we  now  sum  up  the  contributions  of  all  the  modes,  we  obtain  the  final 
solution  [25]  for  the  field  radiated  in  the  shadow  region  by  an  infinitesimal 
magnetic  dipole  of  strength  M located  on  a smooth  convex  conducting  body 
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Q 


Figure  10:  Divergence  of  surface  rays. 


Figure  11:  Geometry  of  the  cylinder  problem. 


2] 


and  the  expression  for  G is  obtained  by  replacing  the  superscript  "h"  by 

"s"  in  (30b),  where  h and  s stand  for  hard  and  soft  boundary  conditions, 

li  s h s 

viz.,  u*0  and  3u/3n  * 0,  respectively.  The  quantities  L ’ , D ’ and 
tl  s 

a ’ , appearing  in  (30b), in  general  depend  upon  the  local  geometry  and  the 
P 

electromagnetic  characteristic  of  the  surface,  frequency  k,  and  the  mode  of 
propagation.  They  are  determined  by  studying  the  asymptotic  expansions  of 
the  exact  solutions  of  some  special  canonical  problems.  Keller  and  Levy 
([20]  and  [27])  have  derived  the  first  few  terms  of  the  asymptotic  expansions 
for  D and  a by  considering  the  canonical  problems  of  scalar  diffraction  by  a 
circular  cylinder,  sphere,  elliptical  and  parabolic  cylinder.  A study  of  the 
above-mentioned  asymptotic  expansions  and  the  works  of  Franz  and  Klante  [28] 
and  Voltmer  [29], who  have  also  investigated  the  same  problem,  as  well  as  a 
comparison  with  the  results  of  the  "direct  integral  equation  approach,"  reveals 
the  following  characteristics  of  the  solution:  i)  the  first-order  terms  in 
the  asymptotic  expansion  of  D and  a are  independent  of  whether  the  problem 


under  consideration  is  scalar  or  vector;  ii)  the  first-order  approximation  of 


D and  a are  dependent  only  on  p^,  the  radius  of  curvature  of  the  surface  along 

the  ray;  iii)  the  second-order  terms  are  functions  not  only  of  p^,  but  also  of 
d2p  d2p 

o , a,  and  o (the  radius  of  the  curvature  of  the  surface  transverse 

~dT~  an 

to  the  ray).  Finally,  the  higher-order  terms  are  different  for  scalar  and 
vector  problems. 


' 
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The  leading  terms  in  the  asymptotic  expansion  of  "diffraction  coefficient" 
D,  "attentuation  constant"  a and  "launching  coefficient"  L are  presented  below: 
"Soft"  polarization: 


. , „1/2  ■ 2‘5/6  • pl/3  • e'3*'12 

[Dp]  * ~YJ 6 ^2 

p k1/6  • [Ai'(-r  )]2 

P 


r • e^^  kp  1/3 

xS  - -B (-=*) 

p P_  2 


L“  - 1 2 (2rk)1/2  A 

P KPa 


2/3 


Ai ' (-r  ) • D 
P P 


"Hard"  polarization: 


[D 


1/2  ,-5/6  A 

h,2  " • 2 • Pa 


1/3 


-jir/12 


p"  k1/6  • r'  • [Ai(-r'  )]2 
P P 


r'  . e^6  kp  , 

•h  - (-^-)1/3 

P P_  2 


(31) 


(32) 


(33) 


(34) 


(35) 


and  Ai(-r  ) - 0,  Ai'(-r')  = 0,  (Ai'  is  the  derivative  of  Ai  with  respect 
P P 

to  its  argument).  Higher-order  terms  in  the  expansion  of  D,  a and  L have 
been  given  in  [25]  and  [30]  and  in  some  of  the  other  works  on  GTD  mentioned 
earlier. 

The  expression  (30)  is  convenient  to  use  in  the  lit  region.  In  the 
shadow  part  of  the  transition  region,  since  the  exponential  decay  of  the  terms 
in  (30)  is  weak.  The  convergence  of  the  series  representation  is  very  slow. 
Furthermore,  the  series  diverges  in  the  lit  part  of  the  transition  region. 
Consequently,  in  these  regions,  it  is  more  reasonable  to  use  an  integral 
representation  for  the  surface  ray  field,  which,  in  our  case,  can  be  expressed 
in  terms  of  Fock  functions  [25]. 

Attempts  have  been  made  to  establish  the  mathematical  validity  of  GTD 
and  to  minimize  its  "nondeductive  parts"  (parts  which  are  based  upon  physical 
intuition  or  the  study  of  the  asymptotic  solution  of  some  simple  problem 
geometrical  concepts  of  different  kinds  of  rays,  diffraction  coefficients, 
attenuation  constants,  etc.).  Kravtsov  [31]  and  Ludwig  [32]  have 
analyzed  the  field  near  the  caustic  surface  (smooth  envelope  of  a family 
of  ray),  and  have  developed  a "uniform  asymptotic  solution"  in  the  sense  that 
it  is  finite  at  the  caustic  and  reduces  to  geometrical  optics  away  from  the 
caustic. 

3.  Generalization  of  GTD  and  Investigation  of  Alternate  Methods 
3.1  Generalization  of  GTD  to  Arbitrary  Surfaces 

Keller's  generalization  of  GTD  for  the  analysis  of  the  field  diffracted 
from  a smooth  convex  object  is  closely  related  to  what  is  known  as  the 
"boundary  layer  technique"  in  the  theory  of  differential  equations  [43].  On 
the  other  hand,  the  "uniform  asymptotic  theory"  is  analogous  to  the  method 
used  by  R.  E.  Langer  and  F.  J.  Oliver  to  find  the  asymptotic  solutions  of  the 
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second-order  differential  equations  near  their  "turning  points,"  which  are 
counterparts  of  the  transition  regions  in  our  case  [33],  l 34 ] , and  [35]. 

The  procedure  is  based  upon  the  generalisation  of  the  geometric  optical 
interpretation  of  the  circular  cylinder  problem.  The  solution  obtained  by 
this  method  involves  some  functions  with  unknown  phase  and  amplitude,  similar 
to  Bessel  and  Hankel  functions.  Since  the  surface  of  a smooth  object  is 
actually  the  caustic  surface  of  diffracted  rays,  the  above-mentioned  formulation 
is  applicable  in  this  case,  too.  Lewis,  et  al.  [36]  have  modified  this 

solution  to  make  it  satisfy  the  boundary  condition  on  a convex  body.  Using 
ray  formalism,  they  have  obtained  an  asymptotic  solution  in  a complicated  form, 
which  they  call  "creeping  wave”  and  satisfies  the  boundary  condition  on  and  is 
uniformly  valid  near  and  away  from  the  surface.  It  should  be  mentioned  that 
the  method  has  been  developed  primarily  for  scalar  diffraction  problems. 

Creeping  waves  that  are  traveling  on  the  surface  of  the  body  generate 
other  kinds  of  diffracted  rays  in  the  presence  of  the  irregularities  in  the 
geometric  or  electromagnetic  characteristics  of  the  surface.  The  effects  of 
discontinuity  in  the  surface  curvature,  its  higher-order  derivatives, 
or  the  surface  impedance  have  been  studied  by  many  authors  [37],  [38], 

[39]  or  [40].  An  exhaustive  study  of  various  diffraction 

mechanisms  and  corresponding  diffraction  coefficients, 

and  constants  associated  with  the  propagation  of  creeping  waves,  has  been 
carried  out  by  Albertsen  [41]. 

At  this  point,  let  us  examine  one  of  the  most  important  features  of  GTD 
and  its  various  modifications.  GTD  formulation  is  essentially  scalar  in  nature 
and  is  heuristic  in  some  parts.  Thus,  when  GTD  is  applied  to  a vector  problem, 
it  is  not  surprising  that  the  coupling  between  various  components  of  the  fields 


are  neglected,  and  each  one  of  them  Is  treated  as  an  uncoupled  scalar  wave. 

The  other  assumptions  in  GTD  are  concerned  with  the  directions  of  these  field 
components  and  the  kind  of  boundary  conditions  they  satisfy  (see  Sec.  2.4). 

As  mentioned  earlier,  non-deductive  parts  of  GTD  are  based  on  asymptotic 
expansions  of  known  solutions  to  some  selected  "canonical'  problems.  Quite 
often  these  canonical  problems  are  not  general  enough  to  fully  and  accurately 
describe  the  local  behavior  of  the  field  for  an  arbitrary  structure.  Finally, 
most  of  the  canonical  problems  investigated  are  two-dimensional  in  nature.  The 
only  exception  to  this  is  the  sphere.  However,  in  so  far  as  the  geometric 
properties  of  the  surface  are  concerned,  the  sphere  is  a very  special  case 
since  its  radius  of  curvature  is  the  same  in  all  directions  and, consequently , 
the  surface  rays  are  torsionless.  Finally,  GTD  fails  when  the  observation 
point  is  located  in  the  transition  regions,  shadow  boundaries  or  in  the  neighbor- 
hood of  a caustic.  In  each  of  these  regions,  one  needs  to  carefully  modify  the 
GTD  formulas  and  often  such  a modification  is  not  too  simple.  Nevertheless, 
in  spite  of  these  difficulties,  GTD  is  recognized  to  be  a powerful  high-frequency 
technique  for  computing  the  leading  terms  of  the  asymptotic  solution.  Two  of 
the  principal  attributes  of  GTD  are  its  simplicity  and  wide  scope  of  application. 

3.2  Spectral  Domain  Approach 

We  now  examine  an  approach  different  from  GTD  which  uses  the  spectrum 

of  the  induced  current,  or  the  expression  for  the  radiated  field,  as  a starting 
point.  In  order  to  gain  a better  insight  into  the  curved-surface  radiation 

and  scattering  problem  and  to  verify  the  basic  assumption  of  GTD,  it  is 

worthwhile  to  consider  such  alternative  approaches,  particularly  if  they  apply 

to  canonical  problems  which  are  more  general  in  nature  than  those  employed  to 

derive  the  GTD  results.  An  example  of  such  a study  would  be  to  consider  the 
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case  of  surface  ray  propagation  with  non-zero  torsion,  a situation  that  occurs 
when  a magnetic  dipole  source  radiates  from  a location  on  the  surface  of  a 
circular  cylinder. 


The  geometry  of  the  problem  is  shown  in  Fig.  11  (p.  21).  The  radius  of  the 
cylinder  is  a and  the  source, which  is  an  infinitesimal  magnetic  dipole  with 
density  M,is  located  at  the  point  Q described  by  the  spherical  polar  coordinates 
(r*a,  9-90®,  ((’••O0).  Each  point  P on  the  surface  of  the  cylinder  is  defined  by  a 
"geodetical  polar  coordinate"  system  (o,B),  where  a is  the  arclength  of  the 
geodesic  connecting  Q to  P and  6 is  the  angle  between  $ (at  point  Q)  and 
geodesic  QP.  The  local  orthonormal  basis  vectors  (a, 6)  are  also  associated  with 
these  two  parameters.  The  observation  point  in  the  far  field  is  specified  by 
its  spherical  polar  coordinates  (r,9,<t>).  The  radiated  field  at  an  arbitrary 
point  can  be  expressed  in  terms  of  two  potentials , $ ani  Y,  which,  in  cylindrical 
coordinates,  can  be  written  as: 
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For  the  problem  under  consideration,  we  can  express  the  spectral  weight 
coefficients  as 
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In  order  to  derive  an  asymptotic  expansion  of  (37)  and  (38),  we  proceed 
as  follows.  As  a first  step,  we  apply  Watsonls  transformation  to  the  infinite 
summation  with  respect  to  n and  employ  appropriate  asymptotic  formulas  for 
Hankel  functions  with  large  order  and  argument  to  derive  the  following  expressions 
for  (37)  and  (38)  under  the  conditions  that  ka  is  large  and  <j>  small  compared  to  ir: 
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where 


fl  = k z + k [p+a(<j>-ir/2)  ] 
z t 

m * (kj.a/2)1^3 


. m($-tr/2) 


fgj  gQ»  g^  ■ Fock's  functions  defined  in  Appendix  A. 


M,  and  M * components  of  M,  (M-n«0) 
z 


Next,  applying  the  "saddle-point"  technique  to  (42)  and  (43)  and  keeping  only 

the  first-order  terms,  the  far  field  can  be  written  in  terms  of  its  components  along 

the  normal  and  tangent  to  the  surface  at  the  "stationary  point"  P as 
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where 

P^:  is  the  stationary  point  of  n which  turns  out  to  be  the 
same  as  the  point  of  diffraction  predicted  by  GTD. 


Cls  - (^)1/3  (W2)  • sin1/3  9 


* radius  of  curvature  of  geodesic  QP^ 


a ■ arc  length  QP^ 

R 3 the  distance  between  the  point  of  diffraction  P^  and 
the  observation  point 


n^  * x &2>  normal  to  the  surface  at  P^ 


The  details  of  the  derivations  of  (44)  and  (45)  are  given  in  Appendix  B. 

Fig.  12  illustrates  the  geometric  meaning  of  some  of  the  parameters 
appearing  in  (44)  and  (45),  for  the  observation  point  is  located  in  the  shadow 
region.  In  this  case,  , which  is  identical  to  £ given  in  (16),  is  the 
reduced  distance  traveled  by  the  surface  ray  before  leaving  the  surface 
tangentially. 

In  the  lit  region,  the  geometric  interpretations  of  a and  £ are  shown 
in  Fig.  13.  The  rays, like  QP^P,  that  do  not  obey  the  generalized  Fermat's 
principle  are  called  "psuedo-rays"  [25].  The  ray  QP^P  appears  to  travel  along 
the  surface  up  to  the  point  P^  and  then  leaves  the  surface  at  P^  tangentially 
in  the  opposite  direction,  to  reach  the  observation  point  P.  It  should  be 
noted  that  formulas  (44)  and  (45)  give  us  the  contribution  of  the  ray  which 
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Figure  12:  Diffraction  of  rays  fay  a cylindrical  body. 


Figure  13:  Diffraction  of  "psuedo-rays . " 


travels  along  the  shortest  path  on  the  surface,  and  thus,  suffers  the  least 
attenuation.  It  is  not  difficult  to  see  that,  in  general,  there  are  infinitely 


r 


many  rays  which  contribute  to  the  total  field  at  any  observation  point.  However, 
their  contributions  are  very  small  compared  to  that  of  QP^  and  their  phases  and 
amplitudes  can  be  determined  in  a similar  manner  (Appendix  B) . Several  other 
remarks  on  the  formulas  given  in  (44)  and  (45)  are  in  order: 

a)  Numerical  results  indicate  that  good  agreement  between  (49), 

(45),  and  the  exact  modal  solution  is  obtained  for  ka  > 10. 

b)  The  zeroth-order  terms  in  the  asymptotic  expansion  of  the  normal 

component  of  the  field  E are  identical  to  those  given  by  GTD; 

-1/3 

however,  the  k terms  derived  from  the  two  approaches  are 
different. 

c)  Tangential  component  of  the  field,  E , given  by  (45),  also 

p 

is  different  from  the  corresponding  expression  based  on  GTD 
by  a multiplicative  polarization  factor.  Specifically, 


(M*  <j>  ) 

■uoTswarap1  (GTD) 


Consequently,  our  results  agree  in  GTD  only  for  the  circumferential 
ray,  i.e.,  for  S = 0. In  addition,  for  an  axial  i.agnetic  dipole 

->■  A A 

(M  • = 0) , GTD  gives  a nonzero  value  for  the  field  in  the  direction 

our  solution  predicts  that  this  field  is  identically  zero,  a result 
which  is  in  complete  agreement  with  the  exact  solution  for  the 
problem. 

In  contrast  to  GTD,  formulas  (44)  and  (45)  are  valid  in 
respect  to  the  location  of  the  observation  point,  be  it  in  the 
lit,  shadow  or  transition  regions.  Although  not  valid  in  the  paraxial 
region  (6*90°),  they  can  be  generalized  to  work  along  this  direction. 
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Finally,  let  us  consider  the  possibility  of  the  generalization  of 
(44)  and  (45)  to  other  convex  surfaces  of  more  general  nature. 

By  "more  general  surfaces,"  we  mean  those  surfaces  which  are  not  substantially 
different  from  cylinders,  some  examples  being  cylinders  with  noncircular 
(elliptical,  hyperbolic  and  parabolic)  cross  sections  and  conical 
surfaces  with  small  apex  angles.  The  key  step  in  a systematic  approach  to 
generalizing  (44)  and  (45)  is  to  use  the  generalized  definition  of  £ given 
in  (16). 

Fig.  14  exhibits  some  initial  results  of  the  generalization  of  these 
formulas  to  the  case  of  a cone.  It  is  evident  that  results  obtained  from 
the  present  approach  agree  quite,  well  with  the  series  solution  which  is  rather 
tedious  and  time-consuming.  We  also  observe  from  Fig.  14c  that  there  is  a 
noticeable  discrepancy  between  the  analytical  solution  and  the  experiment.  Thus, 
within  the  range  of  experimental  error,  our  results  agree  quite  well  with 
those  published  in  the  literature. 

3. 3 Approach  Based  on  an  Asymptotic  Evaluation  of  the  Radiation  Integral 
of  the  Surface  Current 

As  a final  topic,  we  consider  an  approach  based  on  the  asymptotic 
evaluation  of  the  radiation  integral  expressed  in  terms  of  the  induced  surface 
current  which  is  itself  derived  in  an  asymptotic  manner  for  surfaces  with 
large  radius  of  curvature. 

It  was  shown  in  Sec.  2 that  Fock's  theory  can  provide  us  with  an 
expression  for  the  scattered  field  in  the  neighborhood  of  a smooth  convex 
body  illuminated  by  a plane  wave.  Using  this  solution  in  conjunction  with 
the  reciprocity  principle,  we  can  find  the  far  field  radiated  by  a point  source 
located  on  the  surface  of  the  body.  By  generalizing  the  definition  of  £ in 
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Figure  14:  Comparison  between  spectral  domain  results  (UI) , 

modal  approach  (Hughes) , and  experimental  measure- 
ments (Hughes) . The  UI  results  are  derived  from  a 
generalized  version  of  (44)  and  (45)  for  a cone. 
The  Hughes'  results  have  been  reproduced  from  [63] 
and  are  based  on  a modal  series  of  13  terms.  All 
results  are  for  \/2  radial  slot  on  a cone  of 
half-angle  10°. 
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(b)  i - 220°  cut 
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(c)  4>  = 40°  cut  illustrating  the  discrepancy  of 
theoretical  (Hughes)  and  experimental 
(Hughes)  presumably  attributable  to 
experimental  error. 


A 


"1 


Fock's  theory,  we  can  also  write  the  final  result  in  a GTD  format  and  represent 
it  as  a surface  ray.  The  total  field  at  a point  on  the  surface  is  obtained 
by  adding  all  the  possible  rays  which  reach  the  observation  point  P.  Various 
techniques  can  be  used  to  determine  the  field  propagation  along  these  rays. 

For  instance,  when  the  source  is  located  on  the  surface,  and  the  surface  is  a 
conical  one,  the  field  at  each  point  can  be  decomposed  into  two  parts. 

F = F + F2  (47) 

where  F^  is  the  geometrical  optics  field  when  the  observation  point  is  directly 
illuminated  by  the  source,  and  is  the  creeping -wave  contribution  derivable 
via  an  extension  of  Fock's  theory  when  the  point  is  in  the  shadow  region.  The 
other  term,  F2,  is  the  so-called  tip  contribution,  and  can  be  obtained  by 
physical  optics  or  GTD.  Goodrich  et  al.  [42]  have  applied  this  procedure  to 
find  the  radiation  pattern  of  slot  arrays  on  cones. 

The  approximate  induced  surface  current  distribution  can  be  obtained  by 
Fock's  theory,  GTD  [13],  [14],  [16]  and  [25]  or  some  other  appropriate  high 
frequency  technique.  The  induced  surface  current  due  to  a magnetic  dipole  on 
a perfectly  conducting  circular  cylinder  and  cone  has  been  calculated  by  Chang 
et  al.  [44],  and  Chan  et.  al  [45]  whose  procedure  is  based  upon  an  asymptotic 
expansion  of  the  exact  modal  solution  to  the  above-mentioned  problems.  Lee, 
et  al.  [46]  and  [22]  have  treated  the  same  problem  by  a method  based  on 
Fock's  asymptotic  solution  of  the  problem  of  a sphere  [47].  These  expressions 
for  the  current  distribution  can  be  used  in  the  radiation  integral  representation 
of  the  far  field. 

The  numerical  evaluation  of  this  integral  is  a formidable  task,  especially 
when  the  frequency  is  very  high.  Thus,  it  is  highly  desirable  to  have  an 
analytical  and  explicit  formula  for  the  far  field  expressed  in  terms  of  the 
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surface  current.  We  now  discuss  an  approach  for  accomplishing  this  task  and 
examine  the  problem  of  deriving  an  asymptotic  expansion  of  the  far  field 
radiated  due  to  a point  source  located  on  the  surface  of  a smooth,  conducting, 
and  convex  body  of  an  arbitrary  shape. 

Consider  an  arbitrary  smooth  convex  surface  S shown  in  Fig.  15.  Let 
a magnetic  dipole  source  be  located  at  a point  Q on  S.  We  parametrize  the 
surface  S introducing  a "geodetical  polar  coordinate"  system  with  the  pole 
located  at  Q such  that  an  arbitrary  point  on  the  surface  is  defined  by  a 
pair  of  numbers  (a, 8),  where  is  the  arc  length  of  the  geodesic  QP^  and 

is  the  angle  between  QP^  and  some  reference  direction  at  0.  Unit  vectors 
along  the  constant  parameter  curves  o and  8 are  locally  orthogonal.  The 

A A A A 

unit  normal  to  the  surface,  n,  is  given  by  n * o « 8.  An  element  of  length 
in  this  coordinate  system  may  be  written  as 


2 2 2 
ds  * do  + G d6 


The  radiation  integral  for  the  scattered  far  field  can  be  written 


/ 1 id-ii)  • ds 

41T  J ' K- 


where  R is  the  distance  between  any  point  on  the  surface  and  the  observation 


point.  In  the  geodetical  polar  coordinate  system,  we  can  rewrite  a scalar 
component  of  (49),  say  M,  in  terms  of  a double  integral  of  the  following 


general  form 


M ■ / / F(a,S,P) 

D 


/G  da  d8 


where  we  have  assumed  the  following  form  for  the  surface  current: 


J(P1)  - ftP^  exp(-jk^a)  - J^o  + J^S 

■ k - j k2 ; k ■ 2n/X,  k2  > 0,  k2  <<  k 


(51) 


where  f(P^)  is  a relatively  slowly  varying  function  when  k,  the  free-space 

wave  number,  is  large.  This  assumption  is  based  upon  a close  scrutiny  of 

different  asymptotic  formulas  given  for  the  induced  surface  current. 

When  the  observation  point  P is  located  in  the  shadow  reeion,  the 

main  contribution  to  (51)  comes  from  a small  neighborhood  of  the  stationary 

point  of  the  integrand , and  the  stationary  phase  method  for  multiple  integrals 

([48]  and  [49])  is  applicable.  The  asymptotic  expansion  of  (49)  has  been 
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derived  up  to  the  order  k . The  details  of  calculation  have  been  presented 
in  Appendix  C.  The  final  result  is 


E - — 2 (UQ  + Ux  + 0(k“11/6)) 
8tt 

where 
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A,  B,  and  C are  dependent  upon  geometric  properties  of  the  surface  at  the 
stationary  point  which  turns  out  to  be  exactly  the  same  as  the  "point  of 
diffraction"  of  surface  rays.  The  quantities  A,  B,  and  C are  given  by 


where 


. i j.  fa  . _L  + fa  riSS 

A 2G  2 3a  2G 

a a 


+ a60) 2 - (i/2^-2|]  + o (j) 

9 a“ 


B = LSa/G1/2,  C * - 1/p 


p * radius  of  curvature  of  the  geodesic 

a 

Pg  = geodetic  radius  of  curvature 


L36,  L3a  - coefficients  of  the  second  fundamental  form  of  the  surface  (S) 


A geometric  interpretation  of  these  parameters  has  been  illustrated 

Pg  1/2 

in  Fig.  16.  It  is  evident  from  this  figure  that  [R(P+^  ^1  is  simply  the 

g 

divergence  factor  of  the  rays  leaving  the  surface  tangentially  at  the  point 
of  diffraction.  In  using  formula  (56),  we  should  bear  in  mind  that  the 
various  terms  in  and  are  not  of  the  same  order.  For  example,  in  the 


deep  shadow,  is  exponentially  larger  than  Jg. 


The  formulas  given  in  (56)  have  been  tested  and  compared  with  other 
available  solutions.  An  important  conclusion  derived  from  this  comparison 


1 


is  that  although  the  method  of  radiation  integral  is  based  on  less  restrictive 
assumptions,  it  is  perhaps  not  as  useful  as  the  spectral  domain  approach 
because  the  stationary  point  of  the  phase  of  the  integrand  in  (50)  is  of  the 


second  order,  and  hence,  the  asymptotic  expansion  of  this  integral  converges 
rather  slowly  except  when  kp^  is  very  large  ( = 40  or  more). 
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APPENDIX  A:  FOCK  FUNCTIONS 


In  studies  of  radio-wave  propagation  around  the  earth  by  Van  der  Pol, 
Bremmer,  Pryce,  Fock,  and  others,  and  also  the  later  studies  of  diffraction 
of  electromagnetic  waves  by  certain  bodies  of  revolution  ( [51] , [52] , [53] , 

[54] , [55] , [48] , [56] , [57] , [58] , [59] , [60]  and  [50]),  a class  of  universal 
functions  was  introduced  which  can  be  used  to  predict  the  amplitude  and  the 
phase  of  the  reflected  or  diffracted  field  by  smooth  convex  surfaces  [17]. 

An  exhaustive  treatment  of  these  functions  which,  in  general,  are  defined 
as  Fourier  integrals  having  combinations  of  Airy  integrals  in  their  inte- 
grands, has  been  carried  out  by  N.  A.  Logan  [61].  (See  also  Bowman,  et  al. 

[1]  and  Logan  and  Yee  [17]). 

Since  the  first  extensive  application  of  these  functions  to  diffraction 
theory  was  done  by  Fock,  many  authors  named  them  after  him.  Here  we  list 
only  the  most  important  formulas  and  expressions  for  these  functions  without 
going  through  the  details  of  their  derivations.  We  have  followed  Logan's 
set  of  notations  for  these  functions  [61].  However,  since  his  time  dependence 
factor,  exp(-iwt),  is  different  from  one  we  have  used  throughout  this  paper, 
namely  exp(+jwt),  our  expressions,  listed  below,  are  conjugates  of  what 
have  been  presented  in  [ 61] . 

First  we  start  with  general  definitions.  Fock's  most  general  form  of 
the  "Van  der  Pol-Bremmer  diffraction  formula"  is 


V(x,y1,y2,q)  = exp(jir/4) 


|v(t  - y<) 

t 


v'(t)  - qv(t) 
w2(t)  - qw2(t) 


w2(t 


?>> 


dt 


(A.  1) 
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where  w^(t) , w^Ct),  u(t)  and  v(t)  are  Fock.-type  Airy  functions,  defined  as 


u(t)  = /if  Bi(t) 


, v(t)  = /if  Ai(t) 


w1(t)  = u(t)  + jv(t)  , w2(t)  = w^(t) 

We  note  that  and  w2  can  also  be  defined  as  in  Sec.  2.  y>  and 
y<  are  the  larger  and  smaller  of  the  two  numbers  y and  y V(x, y^y^q)  is 
proportional  to  the  attenuation  suffered  by  an  electromagnetic  wave  generated 
by  a source  located  at  reduced  height  y^  above  the  surface  of  a smooth  con- 
vex body,  when  it  reaches  the  observation  point  located  at  reduced  height 
y2  above  the  same  surface,  x is  the  reduced  distance  between  the  source  and 
the  observation  point  along  the  surface,  and  q is  dependent  upon  the  impedance 


of  the  surface.  Let  us  consider  some  useful  limiting  cases. 

When  y^  = y^  = 0,  then  V(x,0,0,q)  is  denoted  by  V^,  where 


VQ(x,q) 


00  e ^xt  w2(t) 

-co  w2(t)  " qw2(t) 
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When  y^  = 0 and  y2  -*■  °°,  then  V -*•  V^(x,q): 


and  also 


vi(x,q)  = ^ ' J_  w2 (t)  - q*2(t)  ‘ dt 


g(x)  = V (x,0)  = — 


1 • f ” e_Jxt 


/tt  ' -°°  2 
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(A. 4) 


(A. 5) 


(A.  6) 


Based  on  Equations  (A. 6)  and  (A. 7),  a class  of  functions  can  be 
defined : 


f(n)W 


!=- ii 
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gc»(I) 


n -jxt  x 

^ • dt  = 

r w2(t)  dxn 


-^.dt.i^I 

w2(t)  dx11 
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where  T is  any  path  in  the  complex  t-plane  which  comes  from  in  a sector 
defined  by  - it  <_  arg(t)  < - -j  and  goes  to  +»  in  the  sector  -- j < arg(t)  < -j. 

In  what  follows,  we  will  give  the  suitable  formulas  for  f(x)  and  g(x)  in 
different  ranges.  Tabulated  values  and  graphs  of  these  functions  can  be 
found  in  [57], [54]  and  [ 61 ] - 

When  x is  very  large  and  negative,  the  following  asymptotic  expansions 
for  f(x)  and  g(x)  can  be  used  [61]: 


f(x)  - + + ... 
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(A. 10) 


(A. 11) 


The  above  formulas  are  valid  and  accurate  for  x <<  -1.  For  moderate 
values  of  x,  namely,  -1  _<  x <_  1,  it  is  difficult  to  find  an  appropriate 
expression.  Although  there  are  some  analytical  techniques  like  "stationary 
phase  method"  or  "Poisson  summation  formula"  which  may  be  used  to  evaluate 
f^  and  g^  for  these  values,  another  possible  way  which  is  probably  easier 
and  more  efficient  is  to  interpolate  the  tabulated  values  of  these  functions 
in  this  range. 
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In  the  vicinity  of  zero  (|x|  =0),  the  Taylor  expansion  can  be  used 


to  calculate  f and  g.  The  coefficients  are  given  by 

(2/ 3) (n-1/ 4)  - 
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where  t(A,u)  is  the  generalized  "tau"  function: 


oo  n. 

T(X,U)  = l ("  1}  -T  , X > 1 


Vn) 


A2(n) 


B0(n) 


B2(n) 


n=0  (n  + u)' 


1 , Ax(n)  = (n  - 1) 


Jr  2 ...  , 26385 

5 j^5n  - 143n  + — 


/ (29  • 32) 


1 , Bj.Cn)  = -7 (n  - 3/2)/48 

(49n2  + 364n  + 39849/16) / (29  • 32) 


When  x is  large,  and  positive,  residue  series  can  be  used  to 


f(n)  and  g(n): 


f(n)(x)  - gj  (2+7n)ir/6 


. (rp)n  exp (rp  • x • e j5ir/6) 
l Ai’(-  r ) 

p-1  P 
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2m 


(A. 12) 


(A. 13) 


(A. 14) 


compute 


(A. 15) 


APPENDIX  B:  DERIVATION  OF  FORMULAS  (44)  and  (45) 


Here,  we  consider  only  the  derivation  of  the  asymptotic  expansion  of 
$ for  a circumferential  magnetic  dipole.  In  this  case,  may  be  written 


as : 


where 


(2tt)  2 


s<V  = l 

n=-a 


dk 


-jk  z 

e 2 • S(kt) 


i . Hn2)(krp) 
J n<J>  n t 

h<2) (kt.) 


Applying  the  Watson  transformation  to  (B.2), 

- j V (0— 7T) 


s(kt)  - j ■ 


ll<2)(k  0) 

V t 

(2) 

C+D  Hv  '(k  a) 
v t 


Sin  vtt 


dv 


where  C and  D are  shown  in  Figure  17.  Or, 


,C2) 


00  cos  v(ir  - y(ktP) 

_ 

sin  vtt  • Hv  ' (k  a) 
v t 


dv 


s(kt)  = j 


Substituting  the  expansion 

cos  v ( tt  - »)  _ • y y _jv(*i+27Tl) 
sin  vtt  " J if1  e 

where  <j>^  * <j>  and  d>2  “ 2tt  — <p,  in  (B.4),  the  result  will  be: 


S(kt) 


2 00 
■I  I 

i»l  1-0 


- H^2)(ktp)  -Jv(*  +2irl) 

' e • dv 

- H^;(kta) 


(B.l) 


(B.2) 


(B.  3) 


(B.4) 


(B.  5) 


(B . 6) 


Each  term  of  the  above  expansion  is  associated  with  a "creeping  wave" 
travelling  in  a counterclockwise  (i  - 1)  or  clockwise  (i  = 2)  direction 
around  the  cylinder.  Following  the  ray  concept,  each  creeping  wave 
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appears  to  be  travelling  along  a specific  surface  ray.  Now,  as  p + 
(far  zone)  for  each  fixed  v,  we  have  [62] 


00 


-j  (k  p-vir/2-ir/4) 

• e C (B. 7) 

On  the  other  hand,  it  can  be  shown  that  the  significant  contribution 

to  S(kt)  comes  from  a small  neighborhood  of  kta.  In  this  neighborhood,  where 

1/3 

k£a  and  v are  large  and  close  to  each  other  (Ik^a  - vj  < |v|  ),  the  Hankel's 

asymptotic  expansion  (B.7)  is  not  valid  any  longer.  In  this  case,  it  is 
necessary  to  expand  Bessel's  functions  in  terms  of  Fock-type,  Airy  functions, 
w^(t)  and  w^(c),  and  their  derivatives  [16]: 


H(2)(x)  - 
v 


mv^T 


fw2(t) 


— |4t  W (t)  + t2  w'(t) 
60m  1 L 


+ . . .) 


(B.8) 


h(2)'cx)  - 


2 r~ 
m /it 


wl(t)  + 

60m 


4t  w^(t)  + (6  - t ) w ^ ( t ) 


+ . 


(B.9) 


where 


fx 

m = 


1/3 


t = V — X (m  is  very  large) 


Inserting  (B.7)  and  the  first-order  terms  of  (B.8)  and  (B.9)  into  (B.6) 
and  (B.l),  we  obtain 


aieM, 

$ . L. 

2 

(2ir; 


2 » 


I l 

i-1  1-0 


f0(Sil) 


where 

m - (kj.a/2)1/3 

0^  = k^z  + kt[p  + a(<J>i  + 2irl  - ir/2)  ] 
= m(<pi  + 2irl  - tt/2) 


2 


m 


(B. 10) 
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y is  the  path  of  integration  in  the  complex  a-plane,  which  is  shown  in 
Figure  18. 

Now  we  deform  the  path  of  integration  into  the  "steepest  descent  path," 
SDP , passing  through  the  saddle  point  of  the  phase  of  the  integrand.  Per- 
forming the  "saddle'point  integration,"  we  can  derive  the  asymptotic 
expansion  of  (B.15)  for  large  kR^.  The  first  order  term  is: 


. Jir/2  . rka 


l l <cos  S.,)' 

i-1  1*0 


f0(5ils) 


(B. 16) 


where  R ^ an<^  £.qs  are  values  of  these  parameters  at  the  stationary  point 
specified  by  a = 8^-,  • 
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Re  (a) 


r r 


Eqn.  (B.16)is  the  creeping-wave  representation  of  the  far  field.  If  the 
cylinder  is  large  (ka  >>  1)  and  |$|  is  not  very  close  to  it,  then  only  the 
first  term  (i  = 0,  1=1)  has  the  most  important  contribution  to  the  total 
infinite  sum,  and  the  other  terms  are  not  significant.  Neglecting  the  other 
terms,  we  obtain  the  result  given  in  (44  and  45).  It  should  be  emphasized 
that  (44)  and  (45)  are  not  valid  when  ! S ! is  close  to  ir/2  (paraxial  region), 
because  in  this  case,  k^_a  is  very  small,  and  (B.7),  (B.8)  and  (3.9)  no  longer 


apply. 


The  other  formulas  can  be  derived  in  a similar  manner. 
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APPENDIX  C:  ASYMPTOTIC  EVALUATION  OF  THE  RADIATION  INTEGRAL 


Consider  the  following  double  integral: 


U(k) 


;(x,y)  • eJk<Kx,y)  dxdy 


(C.l) 


where  g(x,y)  is  rather  slowly  varying,  and  <j>(x,y)  has  a stationary  point 

(x  ,y  ) inside  domain  D.  The  objective  is  to  derive  an  asymptotic  expansion 
s s 

for  (C.l)  when  k is  large. 

Suppose  g and  <j>  have  the  following  forms  around  (x  ,y  ): 


V1 


V1 


g(x,y)  =■  (x  - xg)  (y  - yg)  g-^x.y),  XQ,  uQ,  > 1 


(C.  2) 


d> (x,y ) = (x  ,y  ) + a (x  - x ) [1  + P(x,y)]  + b.  (y  - y )T[1  + Q(x,y)] 
SS  0,U  S U,T  s 

V 

N.  Chako  [48]  has  derived  the  following  asymptotic  series  for  U: 


U(k)  - Bo  • Z A + «2>«1  + s2>  • 


p,q=0 


X0  + P 


y0  + q 


,,,  sq/t 

0“SO>  0 <kV) 


(C.  3) 


where 


. -JL. . Jk*<vys} 


An/5  y./r  (5x) 

^a.o*  (kbo,x} 


oi^  = exp[jir(AQ  + p)  / (25)  ] , =*  exp  [j7r/(25)] 

Sj_  = exp[jir(u0  + q) / (2t ) ] , S2  = exp  ([jir/(2x)  ] 


(Aq  + p)  (26  + ejlT5) 


(yQ  + q)  (2x  + e'’111) 


- 25 


- 2x 


In  order  to  apply  this  procedure  to  the  integrals  of  the  type  *.0.4)  for 

which 


4>(x,y)  = -ft(a,B)  * -(R  + a) 

(C.4) 

r~r 

g(x,y)  = F(a,B,P)  -g- 

(C.  5) 

When  F is  one  of  the  components  of  J(1  - RR) , one  should  first  determine 
the  stationary  point  of  ft,  wherein  its  first-order  derivatives  vanish.  The 
second  step  is  to  compute  the  various  order  derivatives  of  ft,  J,  R,  . ..,at 
this  point,  and  then  insert  them  into  (C.3).  We  just  give  the  main  formulas 
needed  for  these  derivations. 

Suppose  the  surface  of  the  body,  x(a,B),  is  parametrized  by  a geodetical 
polar  coordinate  system.  As  discussed  previously, in  this  system, a is  the 
arc  length  of  the  surface  geodesic  connecting  the  pole  Q to  x(a,B),  and  B is 
the  angle  between  the  geodesic  and  some  fixed  reference  geodesic  at  Q 
(Fig.  15). 

The  element  of  length  in  this  system  is  given  by 

ds2  = da2  + G(a , B)  dB2  (C.6) 
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Let  us  denote  dx(u)/du  by  x^;  then  we  have  the  following  set  of  relations 


68 


-3G/3a  * . 3G/3p  + , t8B  i 

2 a + 2G  x6  L 3 


■+•  -*■  3G/30  -*■  , T 60 

xBa  = xa6  = “2G~  x8  + L X3 


-x-, 


(C  .7) 


(C.  8) 


(C.9) 


and 

where  x = a,  x./*'^  = 6 are  unit  vectors  along  8 51  const,  and  a - const,  curves, 
o A 8 


and 


= n 


x x x. 
o 8 


& 


(C.10) 


is  the  outward  unit  normal  to  the  surface.  Another  quantity  of  interest  is  the 


"geodetical  curvature"  < given  by 

8 


3G/30 


g 2G 

Using  the  above  relations,  we  can  derive  the  following  expressions 
which  holds  true  at  the  stationary  point: 


(C.ll) 


- 1 - i . 1 - 0 

30  0 


ia.  -r 

36 


where  p 


g 


32n 

302 


1/k  , and 
g 


2 

3zn 

3o38 


32Q 

362 


W 


(C. 12) 

(C. 13) 
(C.14) 


33n 

3o3 


33n 

30238 


-LS<V 


(C. 15) 
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3^  _ _G_  3G/3a  32G/3ct2  _ L^_ 

3823o  R2  R 2 Pa 


(C. 16) 


3 J2  33G/38  . 3G  3G  1 . . 6ct  .86.1 
36  * to  • 4G  + L • L + 2 


32G 

3ct38 


(C. 17) 


where  is  the  radius  of  curvature  of  the  geodesic. 

Equations  (C.12)  and  (C.13)  determine  the  location  of  the  stationary 

point.  At  this  point  R * x^,  which,  if  we  introduce  the  ray  concept, 

tells  us  that  the  surface  rays  leave  the  surface  at  the  "point  of  diffraction" 

tangentially.  Equation  (C.14)  indicates  that  the  stationary  point  is  of 

second  order,  so  that  we  need  higher-order  derivities  of  the  phase . L<Ja,  L^a 
88 

and  L are  coefficients  of  the  second  fundamental  form  of  the  surface 
evaluated  at  the  stationary  point.  They  are  defined  as 


ao  -*■  -*■  _ ctS  ■+• 

L = x • x_  , L ax0 
on  3 ct6 


x3  , L 
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X66  * X3 


Using  the  relationships  given  above,  one  can  find  the  expansion  co- 
efficients g,  , , a , b and  A in  (C.3).  Zeroth  and  first-order  terms  in 
skl  mn  mn  pq 

(C.3)  give  us  formulas  (5.6). 

A few  remarks  should  be  made  concerning  the  expansion  presented  in 
(C.3).  First  of  all,  (C.3)  is  a doubly  infinite  series;  therefore,  for  each 
fixed  power  of  k ^ a finite  number  of  terms  should  be  summed  up.  The 
coefficients  of  various  terms  in  these  finite  sums,  namely  A 's  , become 

pq 

very  complicated  when  p and  q are  greater  than  0 or  1.  Another  difficulty  with 
this  series  is  that  when  the  stationary  point  of  the  phase  is  of  an  order 
higher  than  1,  the  difference  between  the  order  of  the  successive  terms 
(when  they  are  ordered  according  to  the  descending  power  of  k)  becomes  very 
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small,  and  consequently  the  infinite  series  converges  very  slowly.  For 
instance,  in  our  problem  where  <5  = 3 and  x = 2 (stationary  point  is  of 
second  order) , sometimes  the  difference  between  the  orders  of  successive 
terms  is  k , which  indicates  the  weak  convergence  (in  an  asymptotic 
sense)  of  the  expansion  in  the  cases  where  the  frequency  is  not  very  large. 
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